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ABSTRACT 

During the FIRE cirrus IFO, the High Spectral Resolution Lidar (HSRL) was operated from a rooftop site 
on the University of Wisconsin- Madison campus. Because the HSRL technique separately measures the molecular 
and cloud particle backscatter components of the lidar return, the optical thickness is determined independent 
of particle backscatter. This is accomplished by comparing the known molecular density distribution to the 
observed decrease in molecular backscatter signal with altitude. The particle to molecular backscatter ratio 
yields calibrated measurements of backscatter cross section that can be plotted to reveal cloud morphology 
without distortion due to attenuation. Changes in cloud particle size shape and phase affect the backscatter to 
extinction ratio ( backscatter-phase function). The HSRL independently measures cloud particle backscatter 
phase function. This paper presents a quantitative analysis of the HSRL cirrus cloud data acquired over an 
-~33 hour period of continuous near-zenith observations. Correlations between small-scale wind structure and 
cirrus cloud morphology have been observed. These correlations can bias the range averaging inherent in wind 
profiling lidars of modest vertical resolution, leading to increased measurement errors at cirrus altitudes. Extended 
periods of low intensity backscatter were noted between more strongly organized cirrus cloud activity. Optical 
thicknesses ranging from 0.01-1 .4, backscatter-phase functions between 0.02-0.065 sr~* , and backscatter cross 
sections spanning 4 orders of magnitude were observed. The altitude relationship between cloud top and bottom 
boundaries and the cloud optical center altitude was dependent on the type of formation observed. Cirrus 
features were observed with characteristic wind drift estimated horizontal sizes of 5 km-400 km. The clouds 
frequently exhibited cellular structure with vertical to horizontal dimension ratios of 1:5-1 :L 


1. Introduction 

Determinations of the optical properties, structure, 
and the vertical and horizontal extent of cirrus clouds 
have broad applications in remote sensing and the at- 
mospheric sciences. Cirrus clouds reflect incoming solar 
radiation and trap outgoing terrestrial radiation; thus, 
the global energy balance depends upon the optical 
and morphological characteristics of these clouds. 
Scattering and absorption by cirrus clouds affect mea- 
surements made by many satellite-bome and ground- 
based remote sensors. Scattering of ambient light by 
the cloud, and thermal emissions from the cloud, can 
increase measurement background noise. Multiple 
scattering processes can adversely affect the divergence 
of optical beams propagating through these clouds. 

Predicting the effects of greenhouse gasses, aerosols, 
albedo changes, and solar fluctuations on climate re- 
quires the development of models that accurately ac- 
count for the highly variable, nonlinear influence of 
clouds on radiative balance. Models of the feedback 
mechanisms between cirrus clouds and earth’s climate 
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can be tested and improved by studying time series of 
cloud formation, maintenance, and dissipation pro- 
cesses. Good models require good initialization data 
and must generate realistic cloud radiative properties. 
Spatial and temporal histories of the optical and mor- 
phological characteristics of real clouds are required to 
fulfill these needs. 

Because of its precise ranging capabilities, spatial 
resolution, and sensitivity, lidar has played an impor- 
tant role in the detection, depiction, and characteriza- 
tion of cirrus clouds (Evans et al. 1966; Platt et al. 
1987; Sassen et al. 1990), and for the verification of 
cloud heights derived from satellite-bome sensor mea- 
surements (Wylie and Menzel 1989). Lidar systems 
that make one measurement at each range can ade- 
quately determine backscatter intensity distributions 
for optically thin clouds; however, they cannot inde- 
pendently determine absolute optical quantities. As 
optical thickness increases, simple plots of lidar back- 
scatter intensity may produce a distorted representation 
of cloud morphology, and can produce serious errors 
in cloud altitude determinations. This is because the 
lidar return signal from any range depends on both the 
backscatter cross section and the optical depth to that 
range. Single-channel lidar systems may not separately 
measure backscatter and extinction. In order to pro- 
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duce calibrated measurements, single-channel lidar re- 
trieval techniques must be constrained by additional 
information or assumptions about the optical prop- 
erties or distribution of the scatterers (Spinhime et al. 
1980; Klett 1981; Weinman 1984; Elorantaand Forrest 
1986). Direct measurements require the determination 
of a signal intensity calibration at each range. 

The High Spectral Resolution Lidar (HSRL) has 
been specifically designed (Shipley et al. 1983; Sroga 
et al. 1983) to produce calibrated measurements of 
aerosol and cloud particle optical depth, extinction 
corrected cloud morphology, and backscatter phase 
function. The HSRL spectrally separates molecular 
backscatter from cloud particle and aerosol backscatter. 
This is possible because thermal agitation Doppler 
broadens the backscattered spectrum from molecules. 
Small aerosols and cloud particles that contribute sig- 
nificantly to the backscatter are more massive than 
molecules and thus exhibit slower Brownian drift ve- 
locities that produce insignificant broadening of the 
scattered spectrum. Using the known distribution of 
molecular scattering cross section to provide a calibra- 
tion reference at every range, extinction is unambig- 
uously determined from the observed range dependent 
decrease in molecular backscatter intensity. The ratio 
of measured aerosol to molecular backscatter intensity 
provides the aerosol backscatter cross section at each 
range. Because aerosol and gaseous absorption are 
negligible at the HSRL wavelength, the backscatter 
phase function can be directly measured. A more de- 
tailed description of the HSRL theory and the defini- 
tions of lidar measured quantities may be found in the 
Appendix. 

HSRL measurements of cirrus cloud optical prop- 
erties were first acquired (Grund 1987) during the 
FIRE Intensive Field Observations (Starr 1987). In 
this paper we discuss the application of the HSRL to 
cirrus cloud measurements and provide a comprehen- 
sive survey of the cirrus cloud optical properties de- 
termined during the case study period. 

2. System characteristics 

Using a multietalon pressure-tuned Fabry-Perot 
spectrometer, the HSRL simultaneously observes the 
lidar return in two channels. The spectrally narrow 
“aerosol channel,” centered on the transmitted wave- 
length, is most sensitive to aerosol scattering and to 
the central region of the Doppler-broadened molecular 
spectrum. With a prominent notch in the center of its 
bandpass, the spectrally wider “molecular channel” 
accepts the entire Doppler-broadened molecular spec- 
trum, while rejecting much of the aerosol scatter. Thus, 
the signal in each channel represents a different linear 
combination of the aerosol and molecular scattering 
contributions to the lidar return. Complete separation 
of the two channel signals requires the determination 
of a 2 X 2 matrix of inversion coefficients. The two 
coefficients, which represent the aerosol backscatter 


transfer function for each channel, are determined by 
diffusely filling the receiver telescope with attenuated 
laser light and observing the response of the two chan- 
nel signals to a spectral scan of the receiver (Grund 
1987). The coefficients that represent the molecular 
backscatter transfer function for each channel are cal- 
culated from the calibration scan convoluted with a 
model of the molecular backscatter spectrum (Yip and 
Nelkin 1964), corrected for temperature and pressure 
at each altitude. Table 1 summarizes the HSRL system 
characteristics at the time the case study data was ac- 
quired. 

Data acquisition and control of the multietalon 
spectrometer and scanning mirror system are accom- 
plished by a DEC 1 1 /73 computer running a multiuser 
operating system. Real time and time-averaged display 
of range square and log corrected data is provided by 
an in A-scope. Data are acquired in 4 second “shots” 
of 32 000 accumulated, range-resolved lidar profiles. 
Each shot is written to magnetic tape for later pro- 
cessing. 

The tapes are subsequently written to a 2.6 gigabyte 
capacity write once optical disk, which facilitates re- 
peated access to individual segments of large datasets. 
Calibration and data analysis are performed on a VAX 
11/751 computer, which has on-line image display and 
graphics capability. 

Table 1 . Summary of HSRL operating parameters during the 
1986 FIRE cirrus IFO. Because of the narrow spectral bandwidth, 
narrow field of view, high repetition rate, and photon counting de- 
tection scheme, the HSRL was capable of measuring cirrus cloud 
optical properties under day or night conditions, while maintaining 
eye-safe operations. 



HSRL Receiver 


Telescope: 

Primary diameter 
Secondary diameter 
Focal length 
F.O.V. (full width) 

.35 m 
.114m 
3.85 m 
320 mR 

Interference filter 

FWHM 

1 nm at 510.6 nm 

Pre-filter eulons: 

Plate diameter 
Etalon spacers 
Combined FWHM 

50 mm 

1.003, .726 mm 
2.5 pm 

High resolution etalon: 

Plate diameter 
Etalon spacer 
Bandwidth (FWHM) 

150 mm 
12.786 mm 
.6 pm 

Signal detection: 

Photon counting 

>10 mHz 

Photomultiplier tubes: 

EMI Gencom 9863 B/ 100 

low afterpulsing 


HSRL Transmitter 



Laser, Wavelength 
Transmitted beam 
diameter 

Transmitted power 
Bandwidth (FWHM) 
Pulse repetition rate 
Pulse length 

CuG 2 , 510.6 nm 
3u mm 

50 mW at 510.6 nm 
.4 pm 
8 kHz 
15 ns 
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3. Data analysis 

In order to facilitate interpretation of the HSRL 
measurements by a varied community of cirrus re- 
searchers, we have chosen two forms for the represen- 
tation of our data for this case study. The first, grey- 
scale imagery of the uninverted “aerosol channel” 
backscatter intensity, has the advantage of higher tem- 
poral resolution, but lacks a calibrated correction for 
extinction with range. The second data representation 
format, as contour plots of calibrated backscatter cross 
section, reveals the true distribution of backscattered 
intensity, but lacks the ability to depict the fine struc- 
tural detail of the cirrus clouds. This resolution limi- 
tation was imposed by the low average power of the 
laser transmitter that necessitated time averaging to 
mitigate the effects of statistical noise in the inverted 
signals. We have chosen to reserve detailed calibrated 
analysis for this paper while presenting the more qual- 
itative grey-scale imagery in the companion lidar in- 
tercomparison paper (Sassen et al. 1990). 

Although the HSRL was primarily designed for the 
measurement of boundary layer aerosol properties, it 
has been successfully adapted to the task of cirrus cloud 
characterization (Grund 1987). Several difficulties arise 
when assessing cirrus cloud optical properties with the 
HSRL: 

a) the signals are reduced by the. additional range 
to the cirrus clouds ( up to 15 km at zenith in midlat- 
itudes); 

b) the molecular scattering intensity is reduced by 
the lower air density at cirrus altitudes; 

c) low temperatures at cirrus altitudes decrea^ the 
Doppler width of the scattered molecular spectrum 
making accurate separation from the unbroadened 
particulate scattering more difficult; 

d) cloud backscatter intensities can be very large 
compared to molecular scattering thus requiring ex- 
treme accuracy ( ‘-~0. 1%) in the determination of the 
inversion coefficients; and 

e) during daylight operations cirrus clouds scatter 
sunlight into the receiver field of view, increasing the 
background-light-induced noise. Aerosol scattering 
cross sections observed in previous operations rarely 
achieved even thin cirrus scattering cross section values. 

To combat these problems, and still maintain rea- 
sonable calibration accuracy, time resolution, and 
measurement linearity, the following data processing 
algorithm has been applied: 

a) A series of inverted molecular and aerosol profiles 
is generated with ^4 s time resolution. Because drifts 
in the system bandpass and receiver-transmitter tuning 
have nonlinear effects on the relative transmission of 
the receiver channels, a new set of separation coeffi- 
cients is chosen for each of these profiles by the follow- 
ing method: 


1 ) A synthetic calibration scan is produced for each 
shot by a linear interpolation in time between 
the adjacent calibration scans. 

2) The spectra] offset between the center of the re- 
ceiver bandpass and the transmitted wavelength 
is estimated from observations of the ratio of the 
aerosol to molecular channel signals under uni- 
form aperture illumination with transmitted laser 
light. This ratio is measured at ~~2 min intervals 
during operations. The ratio is estimated for each 
shot time by an interpolating spline fit to the 
observations. Since the receiver-transmitter may 
exhibit relative drift of either sign, and the band- 
pass characteristics have a slight asymmetry, the 
sense of the tuning offset is assumed to be the 
same as the overall drift observed between the 
preceding and subsequent calibration scans. With 
the tuning drift determined, the inversion coef- 
ficients are calculated from the synthetic calibra- 
tion scan and a model of the molecular scattering 
spectrum that includes the effects of Brillouin 
scattering. 

b) Because the temporal and spatial distribution of 
extinction is highly variable in cirrus clouds, statistical 
fluctuations in lidar signals cannot be reduced by 
straight forward block averaging of profiles. The reason 
for this is evident in Eq. (A I). Because the signal from 
each range depends on an exponential term that varies 
from profile to profile, simple sums will exhibit a range 
dependent bias toward the shots with least attenuation 
(Milton and Woods 1987). To minimize this effect, 
inverted shots are initially summed only to the extent 
that the backscatter profile can be qualitatively dis- 
cerned from noise (typically ~~30 s). Subjective groups 
of these intermediate sums are formed by adding range 
resolved profiles together, as long as the backscatter 
profiles do not appear markedly different. Optical 
properties are determined from each of these grouped 
profiles, then combined in a time-weighted average to 
achieve the reported time resolution. 

c) In calculating the optical properties, the inverted 
molecular signal has not been applied directly within 
the cirrus clouds. Instead, the HSRL molecular channel 
signal has been smoothed according to the following 
algorithm: 

1 ) Regions dominated by Rayleigh scattering are 
determined both above and below the cloud (this is 
possible because first-cut HSRL backscatter cross sec- 
tion profiles clearly indicate regions of enhanced par- 
ticulate scatter, even in the presence of a small cross- 
talk term). 

2) In these regions, a least squares fit is produced 
from the observed separated molecular signal to the 
expected profile for a pure molecular scattering at- 
mosphere calculated from a radiosonde profile of pres- 
sure and temperature. This is acceptable since the ex- 
tinction contribution from particles is small and most 
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of the signal slope in the separated molecular profile is 
associated with the known decrease in density with 
height. 

3 ) The clear-air observed signals are replaced with 
the smooth best fit estimates above and below the 
cloud. 

4) The cloud optical thickness is determined from 
the decrease in the best fit molecular signal across the 
cloud determined in step 3, while accounting for the 
expected decrease in molecular cross section with al- 
titude. 

5 ) The vertical profile of extinction within the cloud 
is calculated from the separated particulate backscatter 
profile using a Bernoulli solution constrained by the 
optical thickness determined in step 4 and the as- 
sumption of a constant backscatter to extinction ratio. 
The Bernoulli technique is employed in this solution 
to reduce the logarithmic range derivative of the lidar 
Eq. (Al) to a first order equation which can be solved 
explicitly for the extinction profile, given the total op- 
tical thickness and an assumed power law relationship 
between backscatter and extinction (Weinman 1988). 

6) The in-cloud molecular backscatter signal is re- 
placed with a smooth estimate calculated from the ex- 
tinction profile determined in step 5 and the known 
altitude distribution of the molecular scattering cross 
section 0 m (R) [see Eq. (A3)]. In this way, noise is re- 
moved from the molecular scattering profile, while the 
spatial distribution of extinction is closely maintained. 
Backscatter cross sections [Eq. (A5)] are calculated us- 
ing this estimate of P m (R). 

7 ) Backscatter phase functions are reported as bulk 
quantities calculated over the entire depth of cloud be- 
cause the Bernoulli solution employed to mitigate the 
effects of noise assumes the backscatter to extinction 
ratio is constant. 

Penetration of clouds, hence cloud-top altitudes, are 
assured in HSRL data by the presence of molecular 
backscatter signal from above the cloud. Molecular 
signal was clearly evident within the reported time res- 
olutions throughout this case study period, suggesting 
confidence in the cloud-top altitudes indicated at the 
level of the minimum plotted backscatter cross section 
contours. Errors in optical thickness and backscatter 
phase function were determined by estimating the ef- 
fects of photon counting statistics on both the range 
resolved signals and the uncorrelated background signal 
according to the methods detailed in Grand ( 1987). 

Cirrus can sometimes form plate-like crystals that 
fall flat-face-down because of aerodynamics. Conse- 
quently, the HSRL was operated 3°-5° from zenith in 
order to minimize the impact of specular backscatter 
on the measurements of backscatter phase function 
and backscatter cross section. Specular backscatter en- 
hancements by a factor of at least 1 3 have been ob- 
served, and larger enhancements seem possible. The 
e~ l full width of the forward scatter diffraction peak 


was — 1.5° (at 1 .06 ) suggesting a minimum crystal 

diameter of ~~50 *im. Fortunately, at our operating 
wavelength, Ci size parameters are fairly large. This 
reduces the angular width of the diffraction peak so 
that the enhanced backscatter drops off rapidly with 
declination from zenith. 

4. Lidar derived optical properties 

Several periods of cirrus exhibiting differing optical 
and morphological characteristics were observed during 
the case study period. We will begin with a discussion 
of optically thin subvisual cirrus, proceed to a quan- 
titative description of a mesoscale uncinus complex, 
continue with a comparison of altocumulus and ice 
cloud scattering properties, and conclude with the 
characterization of the optically thicker cirrus layer. 
Throughout this discussion we will refer to the optical 
properties defined in the Appendix. In addition, the 
reported optical thicknesses ( r ) have had the effects of 
molecular extinction removed so that they represent 
only the attenuation due to particle scattering. Like- 
wise, backscatter cross sections represent aerosol scat- 
tering quantities without molecular scattering contri- 
butions, and will be denoted ft, w /47r. Because we fre- 
quently observe correlations between small-scale wind 
features and lidar backscatter, isotachs are included on 
each backscatter cross section plot for reference. For a 
pictorial overview of the HSRL observations, the reader 
is directed to the lidar intercomparison paper elsewhere 
in this volume (Sassen et al. 1990). 

a. Subvisible and background cirrus 

The occurrence of subvisible cirrus has been a topic 
of interest in the last few years because of the potential 
effects on IR wavelength remote sensors. In addition 
to the attenuation provided by the cloud, thermal 
emissions from the cloud and ambient light scattered 
by the cloud can contribute to measurement back- 
ground noise. Even tenuous clouds can develop sig- 
nificant optical thickness when probed by shallow angle 
or horizontally viewing long range sensors. In addition, 
low-level aerosol backscatter can complicate the re- 
trieval of optical depth from simple lidar observations 
because calibrations often rely upon the determination 
of regions characterized by pure molecular scattering. 
Goud visibility is a complicated property that depends 
on many factors including contrast, discernible struc- 
ture, and sun angle in addition to optical thickness. 
However, iu this paper, we will adopt the convention 
of Sassen et al. ( 1989) and consider clouds with a zenith 
angle r < 0.03 as subvisible. Within this case study 
period, several subvisible clouds were observed by the 
HSRL. 

Figure 1 shows a contour plot of the /W4* of an 
isolated subvisible cirrus cloud. The contour interval 
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F’C. I. Contour plot of the backscatter cross section ( I0~ 4 km" 1 
sr in solid line ) of a subvisible cirrus cloud. The dashed line indicates 
the optical midcloud height (see text). The average optical thickness 
of this cloud between 2245 and 2320 UTC was 0.03 with an average 
bulk backscatter phase function of 0.028 sr-'. Isotachs (dotted line) 
are plotted in m s 1 from the Fort McCoy radiosondes launched at 
2100 and 0000 UTC. 


is 5.0 ( 10) 4 km ~ 1 sr 1 starting from 3.0 ( 10 )"“ km' 1 
sr 1 . The dashed line marks the optical midaltitude of 
the cloud; half the cloud optical depth is accumulated 
to either side of this line. This height was determined 
from the Bernoulli solution for the in-cloud extinction 
profile. The maximum value of 0 aw /4r measured in 
this cloud was 4.1 (10)“ 3 km -1 sr -1 at 2301 UTC. 
Assuming negligible temporal evolution, translation 
with the 25 m s 1 ambient wind suggests a horizontal 
extent of ~60 km for this cloud. When the 2245-2322 
data are averaged together, the mean optical thickness 
for this cloud was determined to be 0.03 ± 0.006. The 
bulk-average backscatter phase function for this cloud 
was found to be 0.028 ± 0.007 sr -1 . Backscatter in- 
tensification near 2300, 2310, and 2320 at ~8 km sug- 
gests this cirrus cloud was composed of —5 km wide 
shallow cells or bands spaced ~20 km apart, imbedded 
in a less dense background cloud. 

On two occasions during the case study, we have 
observed extended periods of enhanced backscatter 
with little apparent structure. Because the HSRL tech- 
nique does not distinguish between backscatter from 
small aerosols (e.g., dust, H 2 S0 4 droplets, etc.) and 
backscatter from cloud particles, we cannot be certain 
that these subvisible background veils are cirrus; how- 
ever, because more dense cirrus are observed at the 
subvisible veil altitudes both before and after these ob- 
servations, these features appear to be at least associated 
with cirrus. Regardless of particle identification, the 
existence of extended regions of elevated scatter at these 
altitudes have similar impacts on radiative balance and 
the operation of remote sensors. 

A persistent region of enhanced backscatter of this 
type is evident in Fig. 2 from 8.5-13.8 km, and a two 


layer formation apparent in Fig. 3 between 7.0 and 
1 3.8 km. Because of the small backscatter signals from 
these clouds, the range resolution has been degraded 
to ~900 m and temporal resolution has been reduced 
to 30-60 minutes. In higher time and spatial resolution 
images of the raw lidar data, these veils seem to be 
rather stable and are not seen to exhibit the cellular 
structure present in the subvisible cirrus of the type 
shown in Fig. 1 . It is clear from these data that average 
backscatter cross sections of < 1 ( 10)~ 4 km' 1 sr - ' can 
exist over extended regions even when structured 
clouds are not detected. Evidently, there must be several 
mechanisms for the generation and maintenance of 
subvisible clouds. 

Because radiosonde humidity measurements do not 
often extend above ~10 km and are frequently un- 
reliable at these altitudes, it is difficult to suggest particle 
composition from the availability of moisture for the 
formation of ice particles in these clouds. However, 
radiosondes launched from Platteville at 0000 and 0300 
show consistently moist (RH with respect to ice >70%) 
regions from 8.8 km-9.5 km and 8.2 km- 10.0 km, 
respectively. The 0900 and 1200 soundings indicated 
moist layers 7.2 km-10.2 km and 6.6 km-10.1 km, 
respectively, both with maxima RH with respect to ice 
in excess of 100%. 

Between 0000 and 0400, the subvisible cloud exhib- 
ited an average optical thickness of 0.01 ± 0.004. From 
0920 to 1 140 the average optical thickness was 0.003 
± 0.002. Wind-drift estimates of the size of the features 



FIG. 2. Background veils of enhanced particle backscatter are fre- 
quently observed at cirrus altitudes. The average optical thickness of 
this layer was —0.01. Backscatter cross section contours (solid line) 
are in units of 1 0 4 km 1 sr ', optical midcloud heights are plotted 
in dashed line, and isotachs are plotted in dotted line in m s' 1 from 
the Fort McCoy radiosondes launched at 0000, 0300 and 0900 UTC 
( the 0600 radiosonde was not launched ) . Note the correlation between 
wind speed and backscatter cross section contour patterns. Correla- 
tions between wind velocity and backscatter distribution (e.g., 1000 
UTC between 9 and 10.5 km) can bias wind profiles acquired from 
future space based Doppler lidars with insufficient vertical resolution 
(see text). 
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FIG. 3. A two-layer region of enhanced backscatter. Backscatter 
cross section contours (solid line) are in units of 10 -3 kin" 1 sr~ l , 
optical midcloud heights are plotted in dashed line, and wind speeds 
arc plotted in dotted line in m s * from the Fort McCoy radiosondes 
launched at 0900 and 1200 UTC. Backscatter cross section values 
are —half of those shown in Fig, 2, and although the wind profiles 
are Quite different, correlations between wind and backscatter are 
still evident. 


depicted in Figs. 2 and 3 suggest that they extend for 
more than 400 km and 180 km, respectively. One im- 
plication of such large-scale optically thin clouds is that 
remote sensors attempting to view horizontally within 
such layers could easily encounter significant optical 
thickness, even though the cloud may not produce a 
visual manifestation. Another implication of such ex- 
tensive cloud blankets is they can alter earth’s radiative 
balance, if present on a global scale, while going un- 
detected by space-borne passive remote sensors. 

Isotachs plotted in Fig. 2 show a general correlation 
between wind speed maxima and backscatter cross 
section contours. Note the 26 and 29 m contours 
closely follow the pattern of backscatter cross section 
distribution. These apparent pattern correlations have 
been observed in several cases, though a consistent set 
of correlation characteristics has not been established. 
It is not clear whether the sheer is part of the cirrus 
generation mechanism or if cirrus formation and winds 
are both responding to the same environmental forcing. 

Because wind-profiling Doppler lidars rely on back- 
scatter from naturally occurring aerosols, correlations 
between wind sheer and the distribution of backscatter 
can bias vertically-averaged wind measurements of in- 
sufficient range resolution. As an example, a wind pro- 
filing lidar observing the 0200 profile between 7.5 and 
14.5 km with 1 km resolution would produce an rms 
wind speed error of ~0.4 m s~‘ due solely to range 
averaging with the maximum error of -1.2 m s' 1 oc- 
curring at 9.8 km. When weighted by the backscatter 
cross section at each range, the 1 km averaged wind 


speed estimates produced rms errors exceeding 0.8 m 
s 1 with a peak error of 2.4 m s' 1 at 8.6 km. 

Figure 3 shows another example of a strong corre- 
lation between wind sheer and backscatter cross section. 
Note the packing of isotachs between 9.5 and 10.5 km 
coincident with the thin cloud region at 1000 UTC. 
When subjected to the above wind profiling lidar anal- 
ysis, this data produced rms errors of —0.5 ms -1 with 
a peak error of 1.7 m s" 1 at 9.8 km due solely to the 
1 km range averaging of the wind speeds between 6.5 
and 14.5 km. Backscatter weighting of the range av- 
erage produced an rms wind speed error of ~1.0 m 
s 1 with a peak error of -4.1 m s -1 at 9.0 km. A study 
of the effects of velocity-backscatter correlations and 
resolution specifications on wind profiling lidar mea- 
surements is reported in Grund et al. (1990). 

b. Mesoscale uncinus complex (MUC) 

Figure 4 presents a contour map of the absolute 
backscatter cross section from a mesoscale uncinus 
complex (MUC) observed between 0500 and 0900 
UTC 28 October. For the basic structural description 
and a pictorial view of the relative backscatter from 
this mesoscale uncinus complex refer to Sassen et al. 

( 1990). A series of uncinus generating cells is evident 
between 9.5 and 1 1.0 km, particularly between 0530 
and 0700. Each of these cells is about 1 50 m thick. 
Cloud translation with the 10 km wind (20 ± 10 m 
s 1 ) would suggest the cells are ~4- 1 2 km across; thus, 
they have a height to width aspect ratio in the range 
~ 1:50-1:180. Averaging times of ~12 minutes were 
chosen for each profile in this cross section so that ex- 



Fig. 4. Ck. - backscatter cross section ( I0~ 5 km' 1 sr 1 , ) 

and optical midcloud height ( ) of a mesoscale uncinus complex 

The average optical thickness of this system be. . . 0600 and 0750 

UTC was 0.58 ± 0.05 which varied from 0.09 + 0.03 at 0750 L fC 
to I . I ± 0.3 at 07 1 8 UTC. The bulk backscatter phase function av- 
eraged over the same time period was 0.042 ± 0.015 sr" 1 The MUC 
passed over Madison just ahead of a mesoscale wind jet. Again note 

the apparent association of wind speed (ms' 1 , ) with citrus 

morphology. Isotachs were determined from the Fort McCoy radio- 
sondes launched at 0300 and 0900 UTC. 
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pected signal-noise-induced errors in the average cloud- 
backscatter cross section were limited to ±15%. The 
maximum j8 flI /4tr for this MUC was determined to be 
0.024 km 1 sr 1 at 0722 near 8.7 km altitude, preceding 
the passage of a wind jet maximum of ~34 m s _l . 
Both the Ft. McCoy and Platteville soundings at 0300 
and 0900 (0600 was not acquired ) indicate the regions 
above 7.4 km were consistently moist (>70% RH with 
respect to ice ) and occasionally reached supersaturation 
as high as ~ 108% with respect to ice. The radiosonde 
profiles also show an abrupt decrease in relative hu- 
midity below 7.4 km suggesting the steep contour gra- 
dient at 0700 at 7.2 km altitude is related to rapid ice 
crystal evaporation in the dry environment beneath 
the complex. Analysis of the profiles averaged over the 
0600-0750 period indicates this system had a mean 
optical thickness of 0.58 ± 0.05, which varied from 
0.09 ± 0.03 at 0750 to 1.1 ±0.3 at 0718. The bulk 
backscatter phase function averaged over the same time 
period was 0.042 ± 0.0 1 5 sr “ 1 . 

Figures 5, 6, and 7 are the 0600, 0700, and 0800 
GOES IR images covering the IFO study area that show 
the MUC passing directly over Madison (stationary 
white-on-black square in south-central Wisconsin). 
The relatively bright (cold) complex is embedded in a 
less intense eastward moving cloud band that extends 
~NW-SE from a more extensive cloud shield covering 
northern Wisconsin. Cloud image tracking indicates 
the band is moving eastward at ~23 ms -1 , close to 
the wind-drift velocity at cirrus altitudes. The complex 
seems to be propagating southeastward as the band 
moves east. Radiosonde data show the 8.5 km wind 
backing from ~300° to — 270° as the complex passes; 
however, the winds show a slight but consistent veering 
with height above 8 km throughout this period. Figure 
4 indicates the most optically dense cloud regions were 
near 8.5 km. Thus, the apparent southeastward prop- 



Fig. 6. The 0700 UTC GOES IR image shows the mesoscale un- 
cinus complex (MUC) passing directly over the lidar site as it prop- 
agates southeastward along the eastward moving cloud band. The 
lidar time height cross section through the cloud band and MUC as 
it passed overhead was shown in Fig. 4. 


agation of the complex seems to be a displacement of 
the MUC generating region along the west wind-driven 
cloud band, rather than a translation of the complex 
with the ambient wind. 

Isotachs, interpreted from the 0300 and 0900 Ft. 
McCoy radiosonde data, show an apparent relationship 
between wind speeds and the MUC backscatter cross 
section distribution. Note the close proximity in time 
and altitude between the wind speed maximum and 
the backscatter cross section maximum. Also note that 
the cloud bottom occurs where the time height cross 
section indicates relatively steady winds. Further, the 
decrease in cloud top altitude seems to be related to 
the development of a wind minimum at 10.5 km. The 
jet and the wind minimum are of smaller scale than 



Fig. 5. The 0600 UTC GOES IR image showing the cirrus cloud 
band over Madison (white pixel in south-central Wisconsin) ex- 
tending NW to SE from a larger cirrus shield covering northern WI. 
The relatively bright cloud regions just to the NW of Madison moved 
rapidly to the SE along the cloud band over the next two hours (see 
Rg. 6). 



Fig. 7. The 0800 UTC GOES IR image shows the apparent dis- 
sipation of the MUC as it continues to propagate to the SE along the 
cloud band. The lidar (see Fig. 4) shows the backscatter structure of 
the cloud band following the MUC. Displacement estimates from 
Figs. 5-7 suggest the cloud band is moving eastward at ~23 ms - ', 
in agreement with rawinsonde-determined winds at the lidar-deter- 
mined cirrus altitudes. 
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the FIRE IFO radiosonde network. They appear 
strongly only in the Ft. McCoy soundings, and are just 
discernible in the Platteville sounding data; thus, no 
attempt has been made to interpolate the radiosonde 
data to Madison observation times. Future observa- 
tions of this type would benefit from a more dense 
temporal and spatial net of wind observations. 

The MUC observations also demonstrate the limi- 
tations inherent in interpreting cirrus cloud morphol- 
ogy strictly from zenith time height cross section mea- 
surements. If the wind sheer is deduced from the ap- 
parent slope of the virga trails, the wind maximum 
would be expected near the upper cloud level at — 10 
km. In fact, the wind contours indicate the maximum 
of -—35 m s _l near 8.5 km, decreasing to ^18ms _l 
at 10 km. Clearly, the time- height cross section of this 
complex does not represent a stationary phenomena 
translating with the ambient winds. Rather, in addition 
to translation, the clouds are undergoing significant 
evolution. Cloud genesis is probably related to circu- 
lations about the mesoscale jet, and the apparent shape 
of virga may be partially governed by the local distri- 
bution of moisture. Clearly, serious study of the mor- 
phology of such mesoscale phenomena can greatly 
benefit from the acquisition of real time three-dimen- 
sional lidar observations. 

c. Altocumulus versus cirrus 

From 1200 to 1500 UTC 28 October, a two-level 
cirrus-altocumulus formation was observed. Figure 8 



1100 1200 1300 1400 

UTC (10/29/66) 

Fig. 8. Time height backscatter cross section ( 10 -3 km' 1 sr~\ 
interval: 7 ( 10 )~ 3 km' 1 sr*‘) of a two-layer cirrus formation with 
altocumulus embedded in the lower layer. The l ( 10 )“ 3 km" 1 sr H 
contour near the ACu outlines what appears to be ice virga. The 
largest backscatter cross section observed during the case study ( 0. 1 8 
km -1 sr _l ) occurred —1330 UTC at —7.6 km. A table of localized 
backscatter phase function and optical thickness illustrating the dif- 
ferences in optical properties of ACu, ice virga, and the upper cloud 
layer may be found in Sassen et al. ( 1990). Isotachs are plotted in 
dotted line inms’ 1 from the Forth McCoy radiosondes launched at 
0900, 1200, and 1500 UTC 



Fig. 9. Time history of the total cloud optical depth { ) and 

bulk normalized backscatter phase function (- * - • , sr" 1 ) coincident 
with the backscatter cross sections plotted in Fig. 8. Ten minute 
averaging has been applied (see section 3 ) to smooth the data. The 
largest values of backscatter phase function values in this record co- 
incide with the cirrus regions containing altocumulus cells. 


highlights the dissimilarity in lidar backscatter between 
these cloud types. The lower layer cloud appears to be 
formed from a series of relatively dense altocumulus 
cells, each with a 0.5- 1.0 km long ice crystal fallout 
tail. The identification of these clouds as ACu is inferred 
from an all-sky photograph taken at Wausau and the 
fact that similar cloud structure was observed at Ft. 
McCoy and Wausau lidar data around this time period 
(Sassen et al. 1990). Ice virga is inferred from the en- 
vironmental temperature (-25° to -30°C) and from 
grey-scale images of lidar backscatter. Virga appears in 
these images as finely detailed striated sheets that de- 
viate from vertical following the altitude profile of hor- 
izontal wind sheer. Translation with the wind speeds 
at 7.7 km would indicate the ACu were about 20 km 
across, yet they were "-200 m thick, producing a height 
to width ratio of only 1:100. As the ACu developed. 
Fig. 9 shows an increase in both the total cloud optical 
thickness and the backscatter phase function, while the 
optical midcloud altitude decreased from — 1 1 km to 
"-7.5 km. Note that the average cloud backscatter 
phase function peaks with the passage of each ACu 
cell. Because similar cloud features were also observed 
to produce little depolarization at Wausau and Ft. 
McCoy, we can infer, with some confidence, that the 
backscatter phase function peaks suggest these clouds 
are composed of supercooled spherical water droplets 
(see Sassen et al. 1 990 for depolarization data and spa- 
tially localized backscatter phase functions). 

The ACu bands at "—7.6 km produce significantly 
greater backscatter than either the 9-12 km cirrus or 
the lower level virga between 6.4 and 7. 1 km. The /3 ar / 
4 7T maximum reached at 1419 exceeded 0.15 km -1 
sr -1 , nearly an order of magnitude greater than the 
largest j3 fl */4ir observed during this case study in any 
cirrus exhibiting significant vertical development. One 
reason for this is that ACu are primarily composed of 
supercooled liquid water droplets; because of optica! 
surface wave phenomena, such droplets are more ef- 
ficient at backscattering than nonspherical ice crystals. 
A second reason is that supercooled water droplets are 
necessarily quite small; whereas, once glaciation begins, 
cloud particles accumulate mass rapidly because of the 
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reduction of saturation vapor pressure with respect to 
ice. Hence, ice particle fall speeds rapidly exceed ver- 
tical atmospheric motion, and particle mass distributes 
over an extended column, decreasing the local back- 
scatter intensity. The backscatter cross section will also 
diminish as growth processes redistribute water mass 
from many small particles to fewer but larger ice crys- 
tals with substantial fall speeds. This is because particle 
mass increases as the cube of the radius while the scat- 
tering cross section approximately follows the radius 
squared. 

The ACu, when present, dominate the optical thick- 
ness profile. This effect is evident in the dramatic de- 
crease of 4 km in the optical midcloud altitude between 
1200 and 1215. Clearly, the occurrence of ACu can 
have significant impact on the radiative properties of 
the cloud independent of cloud top altitude and bottom 
altitudes. 

d. Cirrostratus 

From 1500 UTC 28 October to the end of the case 
study period, the cirrus observed over Madison as- 
sumed a distinct pattern of vertically developed cells 
with regions of enhanced backscatter appearing near 
both cloud top and cloud bottom (see Fig. 10). The 
vertical elongation of backscatter cross section contours 
suggests an upward and/or downward redistribution 
of cloud particles. The top and bottom cellular mor- 
phology of this layer is most evident ~ 1 6 10 where an 
apparent generating region near cloud top produces a 



Fig. 10. From 1500 UTC until operations were terminated, a cirro- 
stratus layer was observed over Madison which contained a distinct 
pattern of vertically developed cells evident in this time height back- 
scatter cross section ( 10" 3 km’ 1 sr H , ) profile. Regions of in- 

tensified backscatter are found near both cloud top and cloud bottom. 

The optical midcloud height ( ) tends to follow these regions 

independent of cloud top and bottom altitudes. The relatively cloud 
free period — 1900 UTC corresponds to the clear arc region noted 

in satellite imagery by Starr and Wylie (1990). Isotachs ( ) 

are plotted in m s‘‘ determined by radiosondes launched from Ft. 
McCoy at 1500, 1800, 2100, and 2400 UTC. 



1500 1600 1700 1800 1900 2000 2100 2200 2300 2400 

UTC 28 Oct 1986 


Fig. 11 . Time history of the total c’ uid optical depth ( ) and 

bulk normalized backscatter phase function ( , sr _l ) coincident 

with the backscatter cross sections plotted in Fig. 10. Peaks in the 
backscatter phase function tend to coincide with apparent generating 
regions while minima tend to occur between cells independent of 
cloud optical thickness suggesting different particle microphysics for 
the two regimes. 


virga trail that extends to, and perhaps seeds, a second 
region of strong backscatter near ~7.5 km. The lower 
regions (at 1610 and 1640) also appear to be precipi- 
tating as evidenced by the vertical elongation of the 
cloud bottom backscatter cross section contours. Note 
that the midcloud altitude follows regions of enhanced 
backscatter very closely and that there is no obvious 
relationship between the appearance of these regions 
and either cloud top or cloud bottom altitudes. Thus, 
the structures that dominate cloud optical properties 
are not necessarily indicated by either the cloud 
boundaries or the physical midcloud altitude. 

With the exception of two short periods, near 1 540 
and 1915, the cirrus optical thickness exceeded 0.1 
forming an apparent cirrostratus layer (see Fig. 11). 
The gap in cloud cover at ~1900 corresponds to the 
passage of the clear arc identified from satellite imagery 
in Starr and Wylie (1990). Note that the peaks in 
backscatter-phase function (~ 0.03 8 sr _l ) tend to co- 
incide with the occurrence of the regions of enhanced 
backscatter cross section independent of total optical 
thickness, while the regions between the cells usually 
produce the minima in the backscatter phase function 
values of ~0.02-0.026 sr -1 . This suggests differences 
in the microphysical properties for the particles asso- 
ciated with the cells compared with those of the sur- 
rounding cloud. One explanation may be that the lower 
regions of enhanced backscatter are areas of cloud par- 
ticle generation that may contain an abundance of 
small supercooled water droplets (temperature at 7.5 
km: 30°C; at 9.5: 45 °C). Although the gen- 

eration mechanism is unclear, an abundance of ex- 
ceedingly small particles could also account for the ob- 
served backscatter-phase function effects (as well as 
low depolarization ratios, see Sassen et al. 1990). This 
issue cannot be resolved from in situ probe data because 
present instruments are incapable of resolving such 
small particles. However, Heymsfield et al. ( 1990) re- 
port observing clouds containing supercooled water 
droplets using aircraft borne in situ probes. In regions 
between the cells, it may be that sufficient time has 
elapsed so that the number density of droplets has been 
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substantially diminished by differential vapor pressure 
effects between ice and water, or the smallest of particles 
may have grown by some coalescence phenomena. 

The modulation of the bulk backscatter-phase func- 
tion also suggests that some of the apparent structure 
observed in the backscatter cross section plots is gen- 
erated by microphysical changes and is not simply pro- 
portional to the distribution of cloud extinction cross 
section. This effect should be considered when deduc- 
ing cloud morphological and optical characteristics 
from simple plots of lidar backscatter. 

The backscatter contours in Fig. 10 indicate a large 
scale pattern of intensification with a period of about 
45 minutes. However, an examination of the higher 
time resolution pictures of the raw laser returns reveals 
the clouds are composed of many repeated finer scale 
features. Figure 12 is an expanded scale view of the 
1800 to 2230 time period that shows a striated 5 min 
structure imbedded in the larger scale organization. 
This pattern is most clearly developed between 2 1 30 
to 2230. Wind drift estimates would suggest that these 
features have a horizontal scale of about 9- 1 1 km; thus, 
the height to width aspect ratio of the elementary 
structures is — 1:5; whereas the larger scale field of en- 
hanced backscatter exhibits an aspect ratio of ~ 1 :25. 


S. Summary 

The HSRL has been successfully adapted to the task 
of cirrus cloud optical property measurement. The 
transmitter configuration used to acquire the case study 
data produce ~-50 mW of output power and achieved 
eye-safe, direct optical depth, and backscatter cross 
section measurements with 10 min averaging times. A 
new continuously pumped, injection seeded, frequency 
doubled Nd:YAG laser transmitter has recently re- 
duced time-averaging constraints by a factor of ~ 10, 
while improving the aerosol-molecular signal separa- 
tion capabilities and wavelength stability of the instru- 
ment (Grund and Eloranta 1 990). 

The cirrus cloud backscatter cross sections, optical 
thickness, and backscatter-phase functions have been 
determined for the 27-28 October 1 986 segment of the 
HSRL FIRE dataset. Features exhibiting backscatter 
cross sections ranging over four orders of magnitude 
have been observed within this 33 hour period. Ex- 
tended episodes of low backscatter cross section 
( ~~ 10 -7 m _l sr" 1 ) have been observed between struc- 
tured cirrus clouds events. 

During this period, cirrus clouds were observed with 
optical thicknesses ranging from 0.01 to 1.4. Bulk-av- 



Fig. 12. Intensity modulated height time display of the log of the R 2 corrected HSRL aerosol channel signal. The image has 10 times the 
temporal resolution necessary to produce the calibrated attenuation corrected cross section plotted in Fig. 10. The image reveals an ~5 
minute structure ( — ^9—1 1 km wide) imbedded in the larger scale organization shown in Fig. 10. 
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erage backscatter-phase functions in the range of 0.02- 
0.065 sr _l were also observed. A variety of cirrus cloud 
structures were observed with apparent sizes ranging 
from 5 to 400 km during this period that exhibited 
aspect ratios of from 1 :5 to 1:1. 

The altitude relationship between cloud top and 
bottom boundaries and the optical center of the cloud 
is influenced by the type of formation observed. ACu 
and uncinus generating regions tend to concentrate at- 
tenuation in regions of less than 200 m thickness that 
dominate the extinction profile, even when the gen- 
erating cell caps an extended column of virga. Virga 
exhibits complicated fine scale structure, often lying in 
interleaved, sheared sheets. In cirrostratus, imbedded 
vertically-developed cells frequently occupy a signifi- 
cant altitude range, and create large spatial inhomo- 
geneities in optical properties. 

Cirrus morphology and generation processes appear 
to be related to the wind field, probably responding to 
secondary circulations. Better temporal and spatial 
resolution in wind measurements in future experiments 
would aid the understanding of cirrus generation and 
dissipation mechanisms. 
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APPENDIX 

HSRL Theory 


The basic principles inherent to all monostatic lidar 
systems are embodied in the lidar equation: 


P(R) = Ec£R 


- 2 




4x 


X exp | 

(after Collis and Russel 1976). 
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Here, P(R) is the power incident upon the receiver 
from range R , E 0 is the energy of the transmitted pulse. 
£ = A r (c/2) where A r is the receiver area and c is the 
speed of light. The subscripts “a” and “m" denote 
aerosol and molecular scattering quantities. The term 
aerosol includes both small aerosol and cirrus cloud 
particles. ft(R) and ft„(R) are the respective scattering 
cross sections per unit volume, P a (7r, R)/4tt and 
P m (x)/4ir are the respective normalized backscatter- 
phase functions, ft(R) is the extinction cross section 
per unit volume, and M(R) is the contribution from 
multiple scattering. The entire expression in the first 
set of brackets is the total normalized backscatter cross 
section per unit volume, ft( R). Both ft( R ) and ft( R ) 
are due to the effects of aerosols and molecules. The 
exponential term accounts for the two-way optical 
thickness between the lidar and the backscattering vol- 
ume, where the optical thickness accumulated between 
range ft and R 2 is represented as: 


r = 





( A2 ) 


Equation (Al) demonstrates that the lidar return 
depends upon both the local value of ft( R) and upon 
r between the lidar and R. Only a single measurement 
of P(R) is provided at each range by single-channel 
lidar systems, leading to ambiguities in the direct eval- 
uation of ft or ft. In fact, there are an infinite family 
of ft( R) profiles for which an increase in optical thick- 
ness with penetration depth is just cancelled by an in- 
crease in backscatter with range, rendering the medium 
invisible to single-channel lidar systems of low sensi- 
tivity (Grund and Eloranta 1987). In optically thicker 
cirrus clouds, the multiple scattering contribution 
M(R) can further complicate matters by effectively in- 
creasing P(R) in a way that depends upon the unknown 
spatial distributions of the angular scattering phase 
function and optical thickness between the lidar and 
the sensed volume (Eloranta and Shipley 1982). 

As mentioned, several techniques have been em- 
ployed to untangle ft from ft. Each method requires 
significant assumptions about the distribution of scat- 
tered and about the nature and profile of the particle 
backscatter to extinction ratio. The multiple scattering 
contribution, which can become large in returns from 
cirrus clouds, may be minimized by a narrow field of 
view (fov) design for the lidar system. This solution is 
often difficult to achieve and is therefore not frequently 
implemented; however, narrow fov is a necessary re- 
quirement for the high resolution spectrometer em- 
ployed in the HSRL (see Table 1 ), thus the uncertain- 
ties caused by the effects of multiple scattering processes 
are greatly reduced in the data acquired with this system 
(Grund 1987). 

Because the HSRL separately measures molecular 
and particulate backscatter, two lidar equations may 
be written that are coupled by the common unknown 
integral of extinction. Assuming M(R) is negligible. 
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the molecular and particulate lidar equations may be 
written: 


P„(R)R 2 = Eo£0„(R)j- ex p -2 f 0,(r')dr' 

07T JO 


(A3) 


and 


P a (R)R 2 = EoZPAR) 


P a ( x, R ) 

4 7T 

expf-2 
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where the normalized molecular backscatter phase 
function (P m (x)/4x) has been replaced with its ana- 
lytic value, 3/8x. 

Because the molecular scattering cross section (0 m ) 
is determined by the profile of atmospheric density cal- 
culated from a convenient radiosonde (or from cli- 
matology), (A3) is completely defined, and may be 
solved explicitly for the extinction ft from the measured 
lidar return P m (R). The ambiguity in separating ex- 
tinction from backscatter in single-channel lidar re- 
trievals is eliminated by the HSRL .technique since the 
system is effectively calibrated at each range with the 
return from the known molecular backscatter target. 

The aerosol backscatter cross section is unambigu- 
ously determined from the ratio of (A4) to (A3): 


PAR) 


P«(x, R) 
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= P m (R) 
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By solving (A3), the extinction cross section may 
be directly determined from the slope of the molecular 
backscatter return and a profile of molecular back- 
scatter cross section calculated from the in situ density: 

at P) - _ir^ ,n (^(/?)/? 2 ) </ln(/S m (/?))l 

2 [ dR dR 

(A6) 

In the absence of aerosol and gaseous absorption, ft, 
= ft ~ Pm\ therefore, the normalized backscatter phase 
function is uniquely determined from (A5) and the 
atmospheric density profile: 

p «< *»•*) ._ g (R) 3 PAR) 

4x Pm( ’ 8*[(3,(R)-0 m (R)]P m (R) ■ 

(A7) 
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